The Globular Cluster Populations of Giant Galaxies: Mosaic 
Imaging of Five Moderate-Luminosity Early- Type Galaxies 



> 



Jonathan R. Hargis and Katherine L. Rhode 
Indiana University, 727 East 3rd Street, Swain West 319, Bloomington, IN 47405, USA 
jhargis@astro . indiana. edu,rhode@astro . indiana.edu 

ABSTRACT 

This paper presents results from wide-field imaging of the globular cluster 
(GC) systems of five intermediate-luminosity {My ~ —21 to —22) early- type 
galaxies. The aim is to accurately quantify the global properties of the GC 
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systems by measuring them out to large radius. We obtained BVR imaging of 
four lenticular galaxies (NGC 5866, NGC 4762, NGC 4754, NGC 3384) and one 
+3 elliptical galaxy (NGC 5813) using the KPNO 4m telescope and MOSAIC imager 

and traced the GC population to projected galactocentric radii ranging from ~ 20 
kpc to 120 kpc. We combine our imaging with Hubble Space Telescope data to 
measure the GC surface density close to the galaxy center. We calculate the total 
number of GCs (Nqc) from the integrated radial profile and find Nqc = 340 ±80 
for NGC 5866, iV G c = 2900 ±400 for NGC 5813, N GC = 270 ±30 for NGC 4762, 
N GC = 115 ± 15 for NGC 4754, and N GC = 120 ± 30 for NGC 3384. The 
measured GC specific frequencies are Sn between 0.6 and 3.6 and T in the range 
0.9 to 4.2. These values are consistent with the mean specific frequencies for 
the galaxies' morphological types found by our survey and other published data. 
Three galaxies (NGC 5866, NGC 5813, NGC 4762) had sufficient numbers of 
GC candidates to investigate color bimodality and color gradients in the GC 
systems. NGC 5813 shows strong evidence (> 3a) for bimodality and a B — R 
color gradient resulting from a more centrally concentrated red (metal-rich) GC 
subpopulation. We find no evidence for statistically significant color gradients in 
the other two galaxies. 

Subject headings: galaxies: elliptical and lenticular, cD - galaxies: formation - 
galaxies: individual (NGC 5866, NGC 5813, NGC 4762, NGC 4754, NGC 3384) 
- galaxies: photometry - galaxies: star clusters: general 
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Introduction 



Common to all types of galaxies - from dwarfs to giant cD galaxies - is a population 
of globular star clusters (GCs). Studies of the Milky Way GC populati on have shown GCs 
to be among the oldest kno wn stellar populations (~ 11 — 13 Gyr; iKrauss fc Chaboyer 
20031: iForbes fc Brideesl 120101 and references therein) and as such thes e objects provide an 



observable record of the chemica l and dynamical history of the Galaxy (ISearle fc Zinnlll978 



Zinnlll993l ; lAshman fc Z epf 19 981). As lumin ous, compact objects (average My ~ —7, median 
half light radii ~ 3 pc; lAshman fc Zepj|l998l ) . GCs are easily observable out to large distances 
(~ 20—30 Mpc) and a re therefore important p robes of galaxy formation and evolution outside 
the Local Group (see lBrodie fc Straderll2006l and references therein). 



In recent years, the GC systems of giant galaxies have been investigat ed using high- 



resolution, hierarchical galaxy fo r mation simulations i n a cosmological context (IBeasley et al. 



2002 



Kravtsov fc GnedirJ 120051 : Moore et all 120061 : iBolev et all 120091 : iGriffen et al. 



2010 



Muratov &: Gnedinl 120101 ). Critical to testing such simulations are accurate measurements 
of the global properties of galaxy GC systems: total numbers of clusters, specific frequen- 
cies, spatial distributions, and color (metallicity) distributions. We have been carrying out 
a wide-field imaging survey of GC systems of giant spiral, elliptical, and lenticular galaxies 
outside the Local Group [d ~ 15 — 20 Mpc). The survey design is presented in Rhode & Zepf 
(2001; hereafter RZ01) and previous results are given in Rhode & Zepf (2001, 2003, 2004; 
hereafter RZ01, RZ03, RZ04), Rhode et al. (2005, 2007, 2010; hereafter R05, R07, R10), 
Hargis et al. (2011; hereafter Hll), and Young, Dowell, & Rhode (2012; hereafter Y12). The 
primary aim of our survey has been to characterize the global properties of GC systems of 
giant galaxies. Such measurements - made over a range of host galaxy luminosity, morphol- 
ogy, and environment - provide fundamental observational constraints to galaxy formation 
models. 

Studies of the GC systems of moderate-luminosity galaxies in particular [My ~ —21 to 
—22) often lack sufficient spatial coverage to provide directly measured ensemble GC system 
properties. Giant galaxy GC systems often extend to 50-100 kpc from the galaxy center (10' 
to 20' at 15 Mpc), and therefore wide- field observations are necessary to directly measure 
these properties and trace GC systems over their full radial extent. A significant number of 
galaxy GC systems, covering a wide range in host galaxy luminosity, have been studied with 
the Hubble Space Telescope (HST; Cote et al. 2004, Jordan et al. 2007, Kundu & Whitmore 
2001a,b). However, the limited spatial coverage of the widest-field HST instrumentation 
provides radial coverage of only ~ 11 kpc for galaxies at Virgo Cluster distances (~ 3f4 x 3f4 
for the Advanced Camera for Surveys Wide Field Camera; ACS/WFC). Therefore, only 
lower luminosity galaxies {My ~ —21 or fainter) will have full radial coverage of their GC 
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systems in studies with HST. In contrast, wide- field imaging studies have typically focused 
on the most luminous (My ~ —22 or brighter) galaxies, e.g., NGC 4472, NGC 4406 (Rhode 
& Zepf 2001, 2004); NGC 1399 (Dirsch et al. 2003); NGC 1407 (Forbes et al. 2011); 
NGC 5128 (Harris et al. 2004); M87 (Tamura et al. 2006). Consequently, studies of the GC 
systems of intermediate-luminosity galaxies with adequate spatial coverage remain rare in 
the literature. 

This paper presents the results of wide-field imaging study of the GC systems of five 
intermediate-luminosity (My ~ —21 to —22) early- type galaxies: one elliptical and four 
lenticular (SO) galaxies. The basic properties of the target galaxies are given in Table [TJ 
We obtained BVR imaging using the Kitt Peak National Observatory (KPNO) Mayall 4m 
telescope and Mosaic camera and, where possible, complemented this data with published 
and archival data from the Hubble Space Telescope (HST). The galaxies have been chosen 
both to increase the number of moderate-luminosity galaxies with well measured global GC 
system properties and to fill in existing gaps in our survey. In particular, this study doubles 
the number of SO galaxies in the survey, resulting in comparable numbers of spiral, elliptical, 
and SOs in the sample. Hargis & Rhode (2013, in preparation) will present a multivariate 
analysis that will include the five galaxies in this study as well as the ~ 25 other galaxies 
analyzed for the survey to date. 

While each target galaxy has had a previous GC system study, this work is the first wide- 
field CCD study of NGC 5866, NGC 5813, NGC 4762, NGC 4754, and NGC 3384. Given 
the sparseness of this galaxy's group environment, the lenticular NGC 5866 (My = —21.1) 
has been characterized as a field galaxy (see Li et al. 2009 and references therein). The GC 
system was previously studied by Cantiello, Blakeslee, & Raimondo (2007; hereafter CBR07) 
who used archival HST ACS data to study the color distribution of the GC system. The ellip- 
tical g alaxy NGC 5813 (M v — 22.3) is the second brightest member of the NGC 5846 galaxy 
group (IMahdavi et al.ll2005l). Pr evious studies of the total number of GCs in NGC 5813 were 



done by lHarris fc van den Berghl ( 19811 ) and iHopp et al.l ( 19951 ). lHarris fc van den Bergh 
( 119811 ) used sing le-filter photograph ic plate data to estimate the specific frequency and total 
number of GCs. IHopp et al.l ( 119951 ) used BVR, small-format CCD imaging (three 4'.3 x 2^7 
pointings) of NGC 5813 to estimate the total number of GCs and study the color distribu- 
tion of the galaxy's GC system. The edge-on lenticular galaxy N GC 4762 ( My = —21.2) is 



located on the eastern edge of the Virgo Cluster and was noted by ISandagd (119611 ) as having 
"the flattest form of any galaxy known" . The lenticular galaxy NGC 4754 (My = —20.7) 
form s a wide but non-interacting pair with NGC 4762 ( jde Vaucouleurs fc de Vaucouleurs 
1964T ). Both galaxies were included in the ACS Virgo Cluster Survey (Cote et al. 2004; 
hereafter ACSVCS). The tot al numbers of GCs a nd col or distributions of NGC 4762 and 
NGC 4754 were studied by iPeng et al.l ( l2006aj Jbl. 120081 ) . A member of the Leo I group 
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of galaxies, the lenticular galaxy NGC 3384 (My = —20.5) forms a close pair (in projec - 
tion) with NGC 3379 bu t the t wo galaxies are non-interacting (Ide Vaucouleurs et al.lll976l ). 
Harris fc van den Berghl ( Il98ll ) estimated the specific frequency and total number of GCs in 
their single-filter photographic plate study of the NGC 3379 field. 

This paper is organized as follows. Section [2] describes our observations and data re- 
duction procedures. Section [3] presents our methods and results for the detection of the GC 
systems. Section @] describes our analysis of the GC systems. Section presents our results 
on the global properties of the GC systems. Section [6] gives a summary of our results. 



2. Observations and Data Reduction 

The target galaxies were imaged using the Mosaic camera on the KPNO Mayall 4m 
telescope. The Mosaic camera consists of eight 2048 x 4096 CCDs separated by a small 
gap (~ 50 pixels). On the Mayall telescope, the imager has a 36' x 36' field of view with 
0.26" pixels. The images of four of the target galaxies (NGC 5866, NGC 5813, NGC 4762, 
NGC 4754) were taken over five nights in May 2010. NGC 3384 was imaged in March 
1999 as part of a study of the GC systems of four E/S0 galaxies (RZ04). For each galaxy, 
we obtained several images in each of three broadband filters (BVR), dithering between 
exposures to facilitate cosmic ray rejection and provide sky coverage across the CCD chip 
gaps. NGC 4762 and NGC 4754 are separated by ~ 11' and were included a single Mosaic 
pointing. The images of NGC 3379 from RZ04 also include NGC 3384 because the galaxies 
have an angular separation of ~ 7'. We used the reduced images of the NGC 3379/3384 field 
from RZ04 for our analysis of NGC 3384. Figures [T0 show the field of view of the Mosaic 
pointings for the target galaxies. 

The observation dates, exposure times, and number of exposures used in our final, 
stacked image for each target are listed in Table [2j During the May 2010 observing run, 
night two was photometric. The remaining nights were non-photometric but generally clear 
with occasional clouds on some nights. In order to postcalibrate the data taken on the 
non-photometric nights, we obtained sh orter single exp osures of each target in each filter 



(exposure times of 600 s) and images of lLandoltl ( 119921 ) standard fields on the photometric 



night. Using these we derived a photometric calibration solution consisting of color coeffi- 
cients and zero points. For the March 1999 observing run, the NGC 3379/3384 data were 
taken under non-photometric conditions and postcalibrated using additional observations as 
detailed in RZ04. 
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The data were reduced using the MSCRED routines in IRAF q The data were first 
processed using ccdproc, applying the overscan and bias level corrections. Flat fielding 
was done using master dome flats and median-smoothed master twilight flats. The indi- 
vidual Mosaic extensions were then combined into a single FITS image using the msczero, 
msccmatch, and mscimage routines. We next measured and subtracted a constant sky back- 
ground level from each image. Where smooth sky background gradients were present we 
used the imsurf it routine to fit and subtract a 2nd-order polynomial surface from the back- 
ground, having first masked the chip gaps, stars, and target galaxies to obtain an accurate 
measure of the background variations. Lastly, the series of images for each filter were scaled 
to a common flux level by measuring ~ 20 bright stars on each image. In each filter, the 
images for a particular target galaxy were aligned and a single stacked image was created 
using mscstack with ccdclip pixel rejection and no zero point offsets or additional scaling. 
The mean background level of the image used as the flux scaling reference frame was added 
back to the stacked image. The range of the mean full-width at half-maximum (FWHM) 
of the image point-spread functions (PSFs) for the final images is as follows: l'.'l to l'/2 for 
NGC 5866, 0"9 to 1'.'2 for NGC 5813, l'/2 to I'M for NGC 4762/4754, and 1"3 to I'M for 
NGC 3384. 



3. Detection of the Globular Cluster System 

3.1. Point Source Detection and Aperture Photometry 

In order to detect the GC systems of the target galaxies, we use the procedures de- 
veloped in our previous wide-field imaging studies. We create a galaxy-subtracted image 
by smoothing the final combined image, subtracting it from the original, and restoring the 
constant sky background level. All subsequent steps in the analysis were performed on these 
galaxy-subtracted images. We detected all sources on the images using the daof ind routine 
in IRAF with 4a to 5a detection thresholds in each of the three filters. Regions of the 
image which contained saturated pixels (due to bleed trails from bright stars) or high noise 
levels (along the frame edges or near the centers of the host galaxy) were masked out. The 
source lists in each filter were matched to proudce a final source list per field. The total 
numbers of detected objects in BVR are 3858 for NGC 5866, 10207 for NGC 5813, 4903 for 
NGC 4762/4754, and 3224 for the NGC 3384/3379 field. 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associ- 
ation of Universities for Research in Astronomy, (AURA), under cooperative agreement with the National 
Science Foundation. 
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Because GCs at the distances of the target galaxies will be point sources in our images, 
we removed extended sources from our list of detected objects. We used a graphical software 
routine to select the sequence of point sources in a plot of object FWHM versus instrumental 
magnitude, accounting for the increasing scatter at faint magnitudes. Figure shows an 
example of the point source selection for NGC 5866. After requiring that objects pass the 
point source criteria in each of the BVR filters, we are left with the following number of 
point sources for the target galaxies: 1551 for NGC 5866, 5709 for NGC 5813, 2167 for 
NGC 4762/4754, and 1728 for the NGC 3384/3379 field. 

Aperture photometry was performed on the point source objects in order to obtain cal- 
ibrated magnitudes. Here we describe the calibration procedures for the May 2010 Mosaic 
data. The photometry and calibration for the NGC 3384/3379 data was performed in a sim- 
ilar manner; the relevant details can be found in RZ04 Section 3.3. We measured the point 
source objects using 1 x FWHM apertures and applied aperture corrections to obtain total 
magnitudes. The aperture corrections were determined using the mean difference between 
the total magnitude and the magnitude within an aperture of 1 x FWHM for 15 — 20 bright 
stars on the images. The aperture corrections ranged from -0.443 to -0.222 with errors of 
0.001 to 0.004. Because the final BVR science images and the photometric calibration data 
were taken on different nights (see Section [2]), we determined bootstrap offsets to the instru- 
mental magnitudes of the photometric night. These offsets were determined using the final 
science images and the series of shorter exposure images taken on the photometric night. 
The bootstrap offsets ranged from -0.406 to -0.007 with errors of 0.002 to 0.008. The final 
calibrated magnitudes were calculated by applying the aperture corrections, atmospheric ex- 
tinction corrections, and bootstrap offsets to the instrumental magnitudes determined using 
a photometric aperture with a radius of 1 x FWHM. Finally, we corrected the calib rated 



magnitudes for Galactic extinction using the reddening maps of ISchlegel et al.l (119981 ). Ta- 
bled] lists Ay for the target galaxies. For the NGC 4762/4754 field, we used the extinction 
corrections for the center of the frame for both galaxies. The difference in extinction across 
the field ranged from only 0.01-0.02 magnitudes, and therefore has a negligible impact on 
our selection of GC candidates. 



3.2. Color Selection 

Our list of GC candidates is determined by selecting objects from our point source list 
that have magnitudes and BVR colors consistent with those of GCs at the distances of the 
target galaxies. We adopt the same selection procedures used in our previous GC system 
studies and here only summarize the steps described in detail in RZ01. We assume a range 
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of GC absolute magnitudes based on observations of the Milky Way GC luminosity function 
(GCLF) and other well-studied nearby galaxies ( — 11 < My < —4; Ashman & Zepf 1998). 
We select as GC candidates those point sources that have B — V and V — R colors consistent 
with the Milky Way GCs but extend the range of allowed metallicities ([Fe/H] = —2.5 to 
0.0). Lastly, we also consider the likelihood of a point source being a GC by examining the 
appearance of objects in HST imaging (see Section 14.41) and their radial distance from the 
galaxy center. 

Figures loTTTUl show the locations of the point sources in the BVR color-color plane. 
Because contamination from unresolved background galaxies is problematic in selecting GC 
candidates from ground-based imaging studies, we illustrate the expected colors of galaxies 
of various morphological types from z = 0.0 to 0.7 (see RZ01 for details). The color selection 
process was performed for each galaxy as follows. First, we removed any point sources 
brighter than the brightest expected GC at the distance of the galaxy (My = —11). Second, 
because the faintest point sources will have large photometric errors, we applied a faint-end 
magnitude cut to keep only those GC candidates with magnitude errors in all filters < 0.05 
mag (color errors < 0.07 mag). This cut ensures we keep only the most well-measured GC 
candidates and, because the number of background galaxies grows at fainter magnitudes, 
this also reduces contamination. Because such a cut could eliminate faint but bona fide GCs, 
we examined the spatial distribution of the objects being cut to ensure a lack of spatial 
concentration of objects around the galaxy center. We note also that this magnitude cut is 
accounted for in our correction for magnitude incompleteness by the fitting of the GCLFs (see 
Section |4T5|) . Third, we selected as GC candidates those remaining point sources which have 
colors (including photometric errors) within 3a of the Milky Way BVR color-color relation, 
expanded to include metallicities from [Fe/H] = —2.5 to 0.0. We examined in detail objects 
near the galaxy center which only barely fail the selection criteria, ensuring we did not miss 
objects which should be included in our final GC candidate lists. Lastly, we removed any 
objects from our GC candidate list which are likely background galaxies as determined from 
our analysis of archival HST images (see Section [4.2.31) . 

For the NGC 4762/4754 field, the color selection procedure required additional steps. 
For NGC 4762, we masked a region around NGC 4754 in order to eliminate GC candidates 
likely associated with this galaxy. We removed any GC candidates within a 4' radius of 
NGC 4754, having estimated the spatial extent of the galaxy's GC system using a radial 
surface density profile of GC candidates (see Section l4.2.ip . For NGC 4754 we applied the 
same masking procedures, removing all GC candidates within a 5'7 radius of NGC 4762 
from the NGC 4754 GC candidate list. Lastly, because the galaxies were imaged together, 
we adopted a single faint-end magnitude cut (based on the photometric errors) for both 
galaxies. We used slightly different bright-end magnitude cuts due to the differing distances 
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of the galaxies. Figures |8] and [9] show the resulting color-color diagrams. The magnitude and 
color cuts have been applied across the entire field; no spatial/radial cuts have been applied 
other than the respective galaxy masks. Therefore, because of the large field of view of the 
Mosaic images (~ 36' x 36'), the final GC candidate lists for both galaxies contained objects 
in common. 

For NGC 3384, the color selection followed the same procedure described for NGC 3379 
in RZ04. The very shallow depth of the B image means the detected objects in all three 
filters are preferentially blue. Thus the overdensity of GC candidates in the color-color plane 
was markedly skewed toward the blue side of the mean Milky Way BVR GC relation. To 
pick up all the objects in the overdensity, we used a la width on the red side of the selection 
box and a 3a width on the blue side of the box. In addition, we masked the GC system of 
NGC 3379 in our determination of the NGC 3384 GC candidates. We eliminated all GC 
candidates with a 6' radius of NGC 3379, adopting a mask size that extends to the radius 
where the final, corrected radial surface density profile is first consistent with zero within 
the errors (RZ04). 

After applying the color and magnitude selection steps, the final number of GC candi- 
dates for the target galaxies are as follows: 290 for NGC 5866, 1300 for NGC 5813, 481 for 
NGC 4762, 451 for NGC 4754, and 181 for NGC 3384. 

4. Analysis of the Globular Cluster System 

4.1. Completeness Testing and Detection Limits 

The point source detection and magnitude completeness limits were measured using a 
series of artificial star tests on the BVR sciences images. We first computed a model PSF for 
each image derived from the best-fit average PSF of the frame. In a single test, we added 600 
artificial stars to the image with magnitudes which vary around 0.2 mag of a specific value. 
We next applied the identical point source detection criteria used to detect each galaxy's GC 
system (see Section 13. ip and computed the fraction of detected artificial stars. We ran the 
tests over a 4-5 magnitude range and thus obtained a measure of the completeness (fraction 
of detected point sources) as a function of magnitude in each filter. Table [3] lists the 90% 
and 50% completeness limits for the galaxies in this study. 
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4.2. Contamination Corrections 

The GC system detection process (see Section [3]) significantly reduces contamination in 
our GC candidate lists from bright Galactic foreground stars, spatially-extended background 
galaxies, and objects with colors inconsistent with those of confirmed GCs and their range of 
metallicities. Even so, our GC candidate lists contain some degree of contamination. Below 
we discuss how we quantify the remaining contamination in the GC candidate samples. 

4-2.1. Contamination Estimates based on the Asymptotic Behavior of the Radial Profile 

The large spatial coverage obtained with the Mosaic camera provides a means of di- 
rectly quantifying the level of contamination without additional pointings and/or control 
fields. Wide-field imaging studies of extragalactic GC systems typically show a centrally- 
concentrated profile when plotting the azimuthally averaged surface density of GC candidates 
as a function of projected radius from the host galaxy center. Assuming the spatial profile 
of the GC system is not larger than the imaging area, the radial profile decreases smoothly 
to a mean background level. Because at large projected radii very few bona fide GC will 
be present (as evidenced by the decreasing radial profile), the population will be dominated 
by contaminating foreground Galactic stars and background galaxies. Measuring the mean 
surface density (over a large area) at large projected radii from the host galaxy center thus 
provides a direct empirical estimate of our sample contamination. The radial surface density 
profiles can then be corrected for contamination by subtracting the measured asymptotic 
surface density of contaminants. 

We constructed radial surface density profiles for the target galaxies and measured the 
asymptotic contamination values as detailed below. In all cases we used circular annuli 
and explored the effects of choices of bins sizes and bin center locations on the shape of 
the radial profile. We accounted for missing area in the bins (due to masked regions or 
bins which extend off the images) by using the "effective area", i.e., calculating the area 
of unmasked pixels. After adopting a contamination estimate, we computed the radially- 
dependent contamination fraction as the ratio of the number of contaminating objects in the 
bin (surface density of contaminants multiplied by the effective area of the bin) to the total 
number of GC candidates in the bin. We describe the initial radial profiles and asymptotic 
contamination estimates for each galaxy below. 

NGC 5866. - We constructed the radial profile using bins with 2'5 spacing. Given 
the relatively small number of GC candidates for this galaxy (N = 290) spread over the 
entire Mosaic frame, our experiments with various radial bin widths indicated that wider 
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bins were necessary to ensure a significant number of GC candidates per bin and produce a 
smooth radial profile. The region inward of ~ 1' was masked due to the poorer quality of 
the B band imaging, so the bin centers ranged from 1' from the galaxy center to 28'5. The 
surface density falls to a nearly constant level in the bins outward of ~ 10', indicating that 
we have observed the full radial extent of the GC system. We computed the asymptotic 
contamination correction as the weighted mean of the surface density in the five radial bins 
from r = 14f8 — 24'7, finding a value of 0.21 ± 0.02 arcmin -2 (assuming Poisson statistics). 

NGC 5813. - We constructed the radial profile using bins with 1' spacing. The region 
inward of ~ 0f4 is masked, so the radial bins ranged from 0f4 from the galaxy center to 26.4'. 
The surface density decreases smoothly to a nearly constant level starting at the r = 12.9' 
bin but shows some small fluctuations between 13' and 18'. We calculated the asymptotic 
surface density as the weighted mean of the surface density in the ten bins between 15'and 
24', finding a value of 0.69 ± 0.04 arcmin -2 (assuming Poisson statistics). 

NGC 4*762. - We constructed the radial profile using bins with 1'.25 spacing ranging 
from 0'5 to 28'. The surface density decreased smoothly from the galaxy center, falling to a 
nearly constant level at the ~ 4f9 bin. To avoid possible contribution from the NGC 4754 
GC system, we measured the contamination using the area outside the GC systems of both 
galaxies. Using a single large radial bin with an inner edge at 15' and extending to 28', we 
find a surface density of contaminants of 0.40 ±0.03 arcmin -2 (assuming Poisson statistics). 

NGC 4^54- - We constructed the radial profile using bins ranging from 0'45 to 26'7 
using 0'.75 wide radial bins. The profile decreased smoothly to a nearly constant level at the 
3fl radial bin. Because the GC candidate list for NGC 4754 differs from NGC 4762 (due 
to the different bright-end magnitude cuts), we cannot use the NGC 4762 contamination 
estimate. We avoid any contribution from the NGC 4762 GC system by using a single large 
radial bin outside the radial extent of both GC systems. Using one bin with an inner edge at 
17' and extending to 28', we find a surface density of contaminants of 0.40 ± 0.03 arcmin -2 
(assuming Poisson statistics), consistent with the contamination estimate for NGC 4762. 

NGC 3384- - We constructed the initial radial profile using bins with 1' spacing ranging 
from 0^59 (inner edge of innermost bin) to 18f59 (outer edge of outermost bin). The profile 
shows a smooth decrease to a nearly constant level at the r = 5'1 radial bin. The surface 
density profile showed a very slight but insignificant increase above the errors in the r = 
8'1 radial bin. We calculated an asymptotic contamination surface density of 0.18 ± 0.02 
(assuming Poisson statistics) using the nine radial bins from lOfl to 18^1, adopting the 
weighted mean of the surface density in these bins as our final value. This value is in good 
agreement with the asymptotic contamination surface density estimate for NGC 3379 of 
0.22 ± 0.02 arcmin" 2 from RZ04. 
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4-2.2. Models of Foreground Galactic Contamination 

A primary source of contamination in our wide-field imaging studies are foreground 
Galactic stars (predominately in the halo) that have magnitudes and colors that pass the 
GC selection criteria. We can estimate the degree of this contamination and compare the 
results to our measured asymptotic contamination levels using Galac tic star coun t mod els. 



We used the Besangon Galactic stellar population synthesis models of iRobin et al.l (120031 ) to 
predict the number counts of Milky Way stars in our observed fields. Briefly, the Besangon 
models produce a four component structural model of the Galaxy (thin disk, thick disk, 
spheroid, bulge) using stellar evolution tracks, initial mass functions (IMFs), age ranges, 
and metallicities for each stellar population. The simulations output a catalog of stars which 
pass the user-defined magnitude range(s), color range(s), and area on the sky. 

For each galaxy, we use the Besangon models estimate the surface density of Galactic 
stars which pass our bright and faint V magnitude cuts, B — V and V — R color cuts, and 
cover the identical area of sky of our Mosaic fields. Because the models employ Monte Carlo 
methods to generate the stellar populations, we run the simulations ten times for each galaxy 
and use the mean and standard deviation of the ten runs to estimate the total number of 
contaminating stars. We show an example simulation and color selection for NGC 5866 in 
FigurefTTJ plotting the BVR colors of simulated stars which fall in the magnitude range of the 
NGC 5866 GC candidates over the area of the Mosaic pointing. Comparing the simulation 
to the NGC 5866 observations (Figure [6]) shows the effect of photometric uncertainties in 
our observations. 

The Besangon estimates of the foreground Galactic contamination for the six tar- 
get galaxies are as follows: 0.15 arcmin" 2 for NGC 5866, 0.29 arcmin" 2 for NGC 5813, 
0.12 arcmin" 2 for NGC 4762, 0.12 arcmin" 2 for NGC 4754, and 0.11 arcmin" 2 for NGC 3384. 
The typical uncertainty (standard deviation) in the surface density was 0.01 arcmin" 2 . For all 
galaxies, the simulated Galactic contamination estimates are lower than the values derived 
from observations (using the asymptotic behavior of the radial profile; see Section I4.2.ip . 
This is consistent with our expectation that the additional source of contamination comes 
from unresolved background galaxies. 



4-2.3. Estimates of Background Galaxy Contamination from HST imaging of Mosaic GC 

Candidates 

Because of the superior spatial resolution provided by HST relative to our ground-based 
imaging, GC candidates that are unresolved in Mosaic data may be resolved as background 
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galaxies when imaged by HST. We can therefore obtain an estimate of the level of background 
galaxy contamination by examining GC candidates in archival HST images. We downloaded 
Wide Field Planetary Camera 2 (WFPC2) and ACS/WFC data for the target galaxies from 
the Hubble Legacy Archive (HLA) §, selecting data which were pipeline-reduced, calibrated, 
and combined to produce a single stacked image in each filter (Level 2 data). Table 0] lists the 
proposal ID, the target galaxy for the observations, the total exposure time of the combined 
images, the HST instrument, and filter for the archival HST data analyzed in this study. 
After locatin g the M osaic GC candidates on the HST images, we followed the method of 



Kundu et al.l ( 119991 ) and performed aperture photometry of the GC candidates using 0.5 
pixel and 3 pixel apertures. The ratio of counts in these apertures (larger aperture divided 
by smaller aperture) gives a measure of the degree of concentration of the light profile. 
Relative to compact GCs, background galaxies will have a slightly larger sp atial extent and 



therefore show a higher count ratio than GCs. Following iKundu et al.l (11999). for the WFPC 
observations we use a threshold ratio of > 8 for objects imaged on the WFPC2 WF chips; 
we had no GC candidates imaged on the WFPC2 PC chip. For the ACS/WFC images, we 
used a threshold ratio of > 12. We examine each flagged GC candidate visually in the HST 
images to verify that they have an extended spatial appearance. As noted in Section 13. 2\ 
we remove any probable background galaxies from our GC candidate lists. Typically only 
1-3 objects were confirmed background galaxies and therefore removing them from the GC 
candidate lists has a negligible impact on our final results. In addition, because the HST 
pointings are centered on the galaxy, removing these objects does not affect the asymptotic 
contamination estimates determined in Section 14.2.11 

Using the number of background galaxies and the area of the HST pointing, we can 
obtain a measure of the surface density of contaminating background galaxies in the Mosaic 
GC candidate lists. For NGC 5866, because the central regions of the Mosaic images were 
heavily masked, only six GC candidates were located on the ACS /WFC images. Of these six, 
none of the GC candidates showed a high count ratio and therefore we obtained no estimate 
of the background galaxy contamination. For NGC 5813, we found that three of 55 Mosaic 
GC candidates are galaxies and we find surface density of 0.36 arcmin -2 . For NGC 4762, we 
found that only one of the 19 GC candidates was a background galaxy from which we find a 
surface density of 0.12 arcmin -2 . For NGC 4754, found that three of the 16 GC candidates 
are background galaxies and we find a surface density of 0.31 arcmin -2 . For NGC 3384, 



2 Based on observations made with the NASA/ESA Hubble Space Telescope, and obtained from the Hubble 
Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the 
Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). 
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we adopted a background galaxy contamination estimate of 0.08 arcmin 2 from the RZ04 
study of NGC 3379. RZ04 analyzed four HST/WFPC2 fields of the NGC 3379/3384 region, 
including a pointing centered on NGC 3384 (see Table 5 in RZ04). Because our selection of 
GC candidates around NGC 3384 is identical to the NGC 3379 study, the HST contamination 
analysis is also identical and provides a contamination estimate for our study of NGC 3384. 

In general we find that the background galaxy contamination estimates are lower than 
our empirical asymptotic contamination estimates. This is consistent with our expectation 
that the remaining contamination comes from foreground Galactic stars. 



4.3. Coverage of the GCLF 

For each target galaxy we constructed the observed GCLF (corrected for contamina- 
tion and magnitude completeness) following the procedure outlined in RZ01. We assign the 
GC candidates to V band magnitude bins of a specified width (bin widths of 0.5 mag for 
NGC 3384; 0.4 mag for NGC 5866, NGC 5813; 0.3 mag for NGC 4762, NGC 4754). We ac- 
count for contamination by using the radially-dependent contamination fractions calculated 
in Section I4.2.1[ applying these fractions as a correction factor. Magnitude incompleteness 
at a given V band magnitude is calculated using the artificial star test results (Section 14. ip 
and the range of B — V and V — R colors of the GC candidates. Following RZ01, we use the 
BVR completeness results to derive a total (convolved) completeness which accounts for the 
individual completeness in all three filters. The final, corrected GCLF is derived by dividing 
the total number of GC candidates in each magnitude bin (corrected for contamination) by 
the convolved completeness fraction at that magnitude. 

To determine the observed coverage of the theoretical GCLF, we assumed a Gaus- 
sian functional form for the theoretical GCLF and peak absolute magnitudes of My 



—7.3 for SO /spiral galaxies and My = —7.4 for elliptical galaxies (lAshman fc ZepJ 11998 



Kundu k, Whitmordl2001af ). For the assumed distances of the target galaxies (see Table UJ, 
the peak apparent V magnitudes of the observed GCLFs are V = 25.1,24.1,23.8,22.8 for 
NGC 5813, NGC 4762, NGC 4754, and NGC 3384, respectively. We adopted dispersions 
for the Gaussian fits of 1.2, 1.3, and 1.4 mag, excluding bins from the fit which had very 
low numbers of GCs or were less than 70% complete. The observed and corrected GCLFs 
with theoretical fits are shown in Figure H2J The mean fractional coverage of the theoreti- 
cal GCLF by the observed GCLF for the target galaxies were as follows: 0.196 ± 0.005 for 
NGC 5813, 0.313 ± 0.002 for NGC 4762, 0.351 ± 0.001 for NGC 4754, and 0.438 ± 0.007 for 
NGC 3384. For the four galaxies, varying the choice of magnitude bin sizes and starting 
bin of the GCLF showed a 6% — 8% change in the fractional coverage. We account for this 
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variation in our error estimates on the total number of GCs in Section 15.31 

For NGC 5866, we adopted the observed peak apparent magnitude of V = 23.46 ±0.06 
and dispersion of a = 1.13 ± 0.05 as measured by CBR07. In their HST ACS study of the 
GC system of NGC 5866, CBR07 obtained nearly full coverage of the theoretical GCLF at 
better than 50% completeness, yielding a robust measure of the peak apparent magnitude 
and dispersion of the GCLF. Assuming our adopted distance to NGC 5866 we would expect 
a peak apparent magnitude of V = 23.6 ± 0.1, which agrees within the errors with the 
empirical measurement of CBR07. We fit the corrected, observed GCLF assuming three 
different V magnitude bin widths (0.3, 0.4, 0.5 mag) and fixed values for the peak apparent 
magnitude and dispersion. The mean fractional coverage of the theoretical GCLF by the 
observed GCLF was 0.409 ± 0.004 for NGC 5866. 



4.4. Analysis of Archival HST Observations 

In order to provide improved spatial coverage of the GC systems of the target galaxies, 
we complement our Mosaic imaging with archival and published HST imaging studies. The 
superior spatial resolution of HST allows us to study the GC population closer to the galaxy 
center than the Mosaic imaging, providing data in a region where ground-based observations 
are often limited. However, because the HST spatial coverage is relatively small, these data 
alone are often insufficient to measure the full radial extent of giant galaxy GC systems 
(RZ03; R10; Hll). In this section we describe our analysis of archival and published HST 
ACS data for NGC 5866, NGC 4762, and NGC 4754. The field of view of the ACS/WFC 
observations are shown in Figures Q] (NGC 5866) and [3] (NGC 4762, NGC 4754). For each 
pointing, the galaxies were offset only slightly from the center of the ACS /WFC field of view; 
the observations provide radial coverage out to ~ 2'.7 from the galaxy centers. 

NGC 5866. - CBR07 studied the GC system of NGC 5866 using deep HST ACS/WFC 
observations in the F435W, F555W, and F625W filters (see Table 0}. CBR07 selected GC 
candidates using size, shape, and color information for their detected sources, ultimately 
deriving a list of 109 GC candidates. T he ACS filter magni tudes were converted to BVR 



magnitudes using the transformations of ISirianni et al.l (120051 ). Comparing the CBR07 color 



selection criteria to our Mosaic color selection procedure showed excellent agreement and 
thus no additional color cuts were made to the CBR07 sample. To account for magnitude 
incompleteness, we apply a faint-end magnitude cut at the ~ 100% completeness level (V = 
23.9) based on the completeness testing results described in CBR07. After applying the 
magnitude cut, we are left with 71 GC candidates which we use in subsequent steps in the 
analysis. 



- 15 - 



Because the expected contamination from unresol ved background gala xies brighter than 



V ~ 24 should be relatively small in HST studies (IKundu et al.l 119991 ) . we assumed no 
contamination when constructing the radial profile and GCLF for the NGC 5866 HST data. 
CBR07 estimated a contamination level in their sample of only 6% of the total number of GC 
candidates over the entire magnitude range. Thus our sub-sample of 71 HST GC candidates, 
selected to be brighter than V = 23.9, should have less than 6% overall contamination. Using 
our GC candidate list, we estimated the fractional coverage of the GCLF by the HST data in 
a manner similar to our Mosaic GCLF analyses. The GCLF was fit using a Gaussian with a 
peak apparent magnitude of V = 23.46 and dispersion of a = 1.1 as determined by CBR07. 
We varied the magnitude bin sizes from 0.4 to 0.6 mag and find a mean fractional coverage 
of 0.69 ± 0.01. Varying the choice of dispersion between a = 1.0 and a = 1.2 yielded a 2% 
change in the fractional coverage. We account for the variation in our error estimates on the 
total number of GCs in Section [5731 For the construction of the radial surface density profile 
(see Section ISTTI) . we adopted circular radial bins with a 0'A width ranging from Of 17 to 2'.3. 

NGC 4762 and NGC 4751 - NGC 4762 and N GC 4754 were included as target galaxies 



in the ACS Virgo Cluster Survey flCote et al.l 120041 1 and imaged in the F475W and F850LP 



filters (Sloan Digital Sky Survey g and i, respectively; see Table @J. Our GC candidate lists 
for NGC 4762 and 4754 were taken from the published catalogs of Jordan et al. 2009 (their 
Table 4), who use magnitude, size, and shape criteria to derive catalogs of GC candidates 
for each ACSVCS galaxy. To account for magnitude incompleteness, we apply a faint-end 
m agnitude cut at the ~ 100% completeness level based on the completeness tests described 



in 



Jordan et al.l (120091 ) . The VCS study considers completeness as a function of object size, 
background surface brightness, and magnitude. We computed the mean background surface 
brightness and mean size for the full VCS GC candidate lists and used g band completeness 
curve (Table 2 from Jordan et al. 2009) which corresponds closest to the mean values to 
determine the magnitude at which the completeness falls below 100%. For both galaxies 
the full VCS GC candidate lists were cut at g = 24.95, leaving a total of 95 and 69 HST 
GC candidates for NGC 4762 and NGC 4754, respectively. No magnitude cut was necessary 
at the brightest end of the GC candidate lists (to reject foreground Galactic stars), as the 
VCS GC candidate selections were already consistent with rejecting objects brighter than 
My < —11 at the adopted distances of the galaxies. 

To account for contamination, we used the VCS control fields to estimate the surface 
density of contaminants in the NGC 4762 and NGC 4754 fields, respectively. The VCS uses 
the control fields to estimate contamination by scaling the control field observations to match 
the observations of each VCS target galaxy, analyzing the fields in a manner identical to the 
the GC candidate analysis, and deriving catalogs of point sources and their probabilities of 
being GCs. After applying our GC candidate selection criteria to the control field catalogs, 
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we find a total of 81 and 89 contaminating objects which, over 17 HST/ ACS pointings, 
yields a surface density of contaminants of 0.42 ± 0.05 arcmin -2 and 0.46 ± 0.05 arcmin -2 
(assuming Poisson errors) for NGC 4762 and NGC 4754, respectively. We estimated the 
fractional coverage of the GCLF by the HST data in a manner similar to our Mosaic GCLF 
analyses, correcting the observed GCLF for contamination by applying radially-dependent 
contamination corrections. For NGC 4762, we adopted circular radial bins with a 0'6 width 
ranging from Of 1 to 2f5. For NGC 4754, we adopted circular, logarithmically-spaced bins 
ranging from 0'09 to 2f0 to better sample the GC surface density near the galaxy center. 
We fit the corrected GCLF with a Gaussian function assuming three different g magnitude 
bin sizes (0.4, 0.5, 0.6) and fixed values for the dispersion and peak apparent magnitude. 
We adopted the VCS values for the GCLF dispersions of a = 1.7, 1.1 and peak apparent 



magn itudes of g = 24.6,24.0 for NGC 4762 and NGC 4754, respectively flVillegas et al. 



20101 ). The mean fractional coverage of the theoretical GCLF by the observed GCLF was 



0.599±0.002 for NGC 4762 and 0.859±0.001 for NGC 4754. Varying the choice of dispersion 
by ±0.1 showed a 1%— 5% change in the fractional coverage, and we account for this variation 
in our error estimate on the total number of GCs in Section 15.31 



5. Global Properties of the Globular Cluster System 

5.1. Radial Distribution of the GC System 

The spatial distributions of the GC systems were investigated by constructing azimuthally- 
averaged radial surface density profiles centered on each host galaxy. For galaxies that had 
additional HST data (NGC 5866, NGC 4754, NGC 4762), the HST and Mosaic radial pro- 
files were constructed independently. For each GC system, we binned the GC candidates in 
concentric circular annuli. For both the Mosaic and HST data, we adopted the same radial 
bins chosen when we used the initial radial profiles to determine the contamination level (see 
Sections 14.2. II and 14.41) . The final radial profiles were corrected for missing area, contamina- 
tion, and magnitude incompleteness as follows. For each circular annulus we accounted for 
the missing area by computing the effective area of the radial bin. We correct the profile 
for contamination by adopting the radially dependent contamination fractions, applying the 
values determined in Section 14.2.11 for the Mosaic data and in Section 14.41 for the HST data. 
Finally, we corrected the radial profiles for magnitude incompleteness by dividing each ra- 
dial bin by the fractional coverage of the observed Mosaic or HST GCLF (see Sections 14.31 
and 14.41 respectively). The surface density was then computed using the effective area and 
the corrected total number of GCs in each bin. The uncertainties in the surface density are 
calculated assuming Poisson errors in the number of GCs and contaminating objects. In 
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Tables |5] to [9] we list the final, corrected radial surface density data for the six galaxies in 
this study: mean radius of the bin, corrected surface density of GCs and uncertainty, frac- 
tional area coverage of the annulus, and the data source (either HST or Mosaic). Figures IT31 
to [T7] show the radial surface density of GC candidates as a function of the projected radial 
distance. 

The radial surface density profiles follow the trends observed in our previous wide-field 
imaging studies: a strong overdensity of GCs centered on the host galaxy which decreases 
smoothly as a function of projected radius away from the galaxy. As noted in Section 14.2. H 
the profiles decrease to a constant surface density (consistent with zero after applying the 
corrections), which indicates that we have observed the full radial extent of the galaxy's GC 
system. For the three galaxies with both HST and Mosaic data (NGC 5866, NGC 4762, 
NGC 4754) we find excellent agreement between the radial profiles in the regions of over- 
lapping spatial coverage. In all three cases, the radial profiles continue to decrease outside 
the HST field of view as evidenced by the Mosiac data. To further describe the radial 
distribution of GCs, the corrected radial profiles were fit with both a power-law profile 
(logcTcc — a o + a i logr) and a de Vaucouleurs profile (logocc — a o + air 1 / 4 ). The resulting 
fit parameters are listed in Table [TO The integration of the best-fit radial profiles to estimate 
the total number of GCs is described in Section 15.31 

We estimate the radial extent of the galaxy's GC system by examining the point in the 
final corrected radial profile where the surface density becomes consistent with zero. For 
our survey, we have defined the radial extent of the GC system as the radius at which the 
surface density of GCs falls to zero within the errors and stays consistent with zero (RZ07). 
Although this definition is sensitive to the survey methods (corrections for magnitude in- 
completeness and contamination in particular), it allows for a systematic and internally 
consistent comparison of the projected size of the GC system for galaxies in our survey. For 
the galaxies in this study, the large spatial coverage of the Mosaic observations resulted in 
stochastic variations in the surface density around zero in the outer bins of the radial profile. 
For example, the profile of NGC 3384 shows a small variation above the background in a 
radial bin more than 12' from the main spatial overdensity of the GC system. Accounting 
for the observed variations, we adopt the following radial bins as the radial extent of the 
GC systems: r = 9'.8 (44 ± 11 kpc) for NGC 5866, r = 12f9 (120 ± 14 kpc) for NGC 5813, 
r = 4(9 (27 ±7 kpc) for NGC 4762, r = 3(1 (15 ±4) for NGC 4754, and r = 5(1 (17 ±4 kpc) 
for NGC 3384. Our uncertainty in the radial extent includes the uncertainty in the distance 
modulus and a one bin-width uncertainty in our determination of where the surface density 
of GCs falls to zero. 
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5.2. Color Distributions and Color Gradients 



An important goal of our wide-field imaging survey is to quantify the global color distri- 
bution of the GC systems of giant gal axies. For stellar populations with ages > 2 Gyr, broad- 
band optical colors trace metallicity ( IWortheylll994j ; iBruzual & Charlotll2003l ) such that old, 
metal-poor populations have photometrically bluer optical colors than old, metal-rich popula- 
tions. The presence of broad an d/or bimod al color distributions of GC s ystems in giant galax- 
ies - including the Milky Way ( 



Zinn 



1985) and many giant ellipticals ( iKundu fc Whitmore 
2001al ; Kundu fc Zepj 120071 ; IStrader et all 120071 ) - has therefore been largely interpreted as 
tracing an underlying bimodal metallicity distribution in GC populations (see Brodie & 
Strader 2006 and references therein). Because metallicity differences in the subpopulations 



points to different episodes of star formation, the bimodal color distri 
have been a key differentiator between galaxy formation scenarios ( 


Dutions of GC systems 


Ashman & Zepl 


1992; 


Forbes et al. 


1997; 


Cote et al. 


1998; 


Beaslev et al. 


2002; 


Muratov & Gnedin 


2010 


). Quanti- 



fying the properties of GC system color distributions for a large number of galaxies over a 
range of mass, morphology, and environment provides a key means of testing galaxy forma- 
tion scenarios. 

To investigate the color distributions of the galaxies in our sample, we construct and 
analyze a sample of GCs that have 90% magnitude completeness in all three filters. We 
minimize the effects of contamination by performing a radial cut on the 90% complete 
sample, limiting our analysis to only those GCs that are within the radial extent of the 
GC system. For each galaxy, we analyze the result ing sample for color bimodality using the 
KMM mixture modeling code ( lAshman et al.lll994l ) to test whether the distribution is better 
described statistically with one more Gaussian distributions. We explore both homoscedastic 
(same dispersion) and heteroscedastic (differing dispersion) models. We plot the B — R colors 
for the radially-cut, 90% GC sample as a function of projected radius from the galaxy center 
to explore the possibility of color gradients in the galaxy's GC. Linear least-squares fits to 
these data are used to evaluate the presence of color gradients. 

We discuss the results for individual galaxies here (The results regarding color bimodal- 
ity and color gradients are summarized in Table [TT]) : 



NGC 5866. - For NGC 5866, the 90% complete, radially-cut (r < 11') sample resulted 
in a total of 51 objects, just above the 50 objects necessary for reliable KMM results. The 
color distribution, shown in Figure [TSJ is quite broad and has peaks at B — R ~ 1.0 and 
B — R ~ 1.4. The KMM results indicate that the distribution is likely bimodal at the 
97% (> 2a) confidence level (homoscedastic case), finding a blue (lower metallicity) peak 
at B — R — 1.08 and a red peak (higher metallicity) at B — R = 1.43. The KMM analysis 
finds a split in the GC colors at B — R ~ 1.3, with 76% (N = 39) GCs comprising the blue 
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subpopulation and 24% (N = 12) GCs comprising the red subpopulation. Figure [18] shows 
the KMM results with the 90% complete, radial cut GC color distribution. Separating the GC 
sample at the empirically-determined color separation between blue and red subpopulations 
for elliptical galaxies of B - R = 1.23 (RZ01, RZ04, Kundu & Whitmore 2001a), we find 
similar results: a blue fraction of 71% (N = 36) and a red fraction of 29% (N = 15). CBR07 
also performed a KMM analysis of their B — R distribution of 109 GC candidates. They find 
a likelihood that the distribution is bimodal at the 96% confidence level, with distribution 
peaks at B — R — 1.12 and B — R = 1.35 and an equal number of red and blue GCs in 
each subpopulation. For our analysis of the mass-normalized number of blue GCs (Tbi ue ), 
we average the KMM and empirical split results and adopt a blue GC fraction of 74%. An 
examination of the B — R colors for the radial cut, 90% sample as a function of projected 
radius (see Figure HI?]) showed no statistically significant color gradient within r < 11' (the 
radial extent of the GC system). CBR07 find a B — R color gradient in their sample of HST 
GC candidates within r < 2' resulting from the higher central concentration of the red GC 
population in the inner ~ 1'. 

NGC 5813. - The 90% complete, radial cut sample (r < 13:4) for NGC 5813 contained 
528 objects. Figure I2"U1 shows the B — R color distribution for this sample; the distribution is 
broad and has peaks at B — R ~ 1.1 and B — R ~ 1.35, including a red "tail" of objects with 
B — R > 1.5. The KMM results indicate that unimodality can be ruled out at better than 
the 99.99% confidence level in both the homoscedastic and heteroscedastic cases. For both 
mixture models, the blue peak is found at roughly the same location (B — R ~ 1.1 — 1.2), 
while the location of the red peak differs depending on the dispersion model assumptions 
(B — R ~ 1.5 for the homoscedastic case versus B — R ~ 1.3 for the heteroscedastic case). 
The resulting number of red and blue GCs comprising each subpopulation also changes; we 
find a blue fraction of 83% (N = 440) in the homoscedastic case compared to 63% (N = 331) 
in the heteroscedastic case. Splitting the distribution at the empirically-determined color 
separation of B — R = 1.23 we find a blue fraction of 58% (N = 304), in good agreement 
with the heteroscedastic KMM results. We adopt the mean of the three blue fraction values 
(two KMM models and the empirical split) of 68% as our final estimate of the blue fraction. 

With a large number of GC candidates in NGC 5813 and clear evidence for bimodal- 
ity, we can extend this analysis to look for differences in the spatial distributions of the 
red and blue GC subpopulations. We used the empirically-determined color separation of 
B — R = 1.23 to split the GC candidates in the full 90% sample (no radial cut; iV = 809 
objects) into red and blue subpopulations. Figure I2T1 shows the spatial positions of the GCs 
for each subpopulation with respect to the host galaxy center. The red subpopulation is not- 
icably more centrally concentrated compared to the blue subpopulation. This is consistent 
with similar studies of the spatial distributions of GC subpopulations in other giant galaxies 
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(see Brodie & Strader 2006 and references therein). We constructed radial surface density 
profiles (corrected for contamination, magnitude incompleteness, and missing area) for the 
two subpopulations (bottom panel of Figure |2"T1). The radial profile for the red GC subpopu- 
lation is steeper than that for the blue GC subpopulation and becomes consistent with zero 
(within the errors) at ~ 8' (~ 75 kpc). The blue GC subpopulation, however, is more spa- 
tially extended and the profile does not becomes consistent with zero until ~ 13' (~ 120 kpc) 
- the radial extent of the combined GC profile as determined in Section 15.11 An examination 
of the B — R colors as a function of projected radius (see Figure HH]) shows a statistically 
significant color gradient (slope = A(B — R)/Ar = —0.018 ± 0.006) in the inner 6' (~ 55 
kpc) of the GC system. This results from the decreasing fraction of red GCs as a function of 
radius from the host galaxy center within this region of the GC system. Using the GC B — R 
color-[Fe/H] relation from Rhode & Zepf (2001, 2004; derived from the Milky Way GC data 
of Harris 1996), we translated the GC colors to [Fe/H] and fit the metallicities as function 
of radius. This yields a best-fit metallicity gradient of A[Fe/H]/Alogr = —0.36 ± 0.10. If 
we consider the entire radial extent of the GC system (r < 13'4), we find no evidence of a 
statistically significant color gradient. 

NGC 4^62. - For NGC 4762, we constructed a 90% complete, radial cut sample of 
N = 50 objects. Because the radial extent of the GC system (outer edge at r = 5'.5) 
included fewer than 50 GC candidates (only N = 44), we extended the radial cut to r = 5'7 
to create a sample of exactly N = 50 objects. The B — R color distribution is quite broad 
and shows peaks at B — R ~ 1.25 and B — R ~ 1.4. The radial distribution of these GCs 
(see Figure [T9~j) shows that the B — R ~ 1.4 peak is due to a cluster of 10 GC candidates 
at projected radii r < 2' with 1.35 < B — R < 1.45. These objects are likely a population 
of "faint fuzzies" (FFs): star clusters that are typically old, metal-rich, have slightly larger 
sizes than typical GCs, and are often spatially associated with the host galaxy disk (Brodie 
& Larsen 2002; Peng et al. 2006b and references therein). The KMM analysis of the 
iV = 50 sample shows that we cannot rule out unimodality at a statistically significant level; 
the distribution is only consistent with bimodality at the ~ 54% significance level (mean 
significance level of the heteroscedastic an d homoscedastic m odels). The ACSVCS study of 



the g — z color distribution of NGC 4762 f lPeng et al.ll2006ai ) also did not detect bimodality 



at a statistically significant level in their sample of N = 123 GC candidates. Given the 
lack of strong evidence for color bimodality, we estimate the blue fraction of GCs using the 
empirically determined split at B — R = 1.25. For the N = 50 sample (with radial cut at 
r = 5'.7) we find a fraction of blue GCs of 34% (N = 17). If we remove the ten FFs from 
the sample, we find a blue GC fraction of 43% (N = 17) using the empirically-determined 
split. We therefore adopt a blue GC fraction of ~ 40%. We do not require bimodality to 
estimate a meaningful blue fraction; simply using the empirical split allows us to compare 
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the blue GC population in a self-consistent manner. Examining the B — R colors of our 90% 
complete, radial cut sample as a function of projected radius (see Figure [T9|) . we find no 
evidence of a statistically significant color gradient within r = 5'7 in NGC 4762 (either with 
or without the FF population). The ACSVCS study of g — z color gra dients (limited to the 



inner 2'A) also did not find evidence of a color gradient in NGC 4762 (ILiu et all 120111 ). 



NGC 4754- - For NGC 4754, the 90% sample based on a radial extent of the GC system 
(r = 3'.5) contained only iV = 24 objects. To create a sample of the minimum N = 50 objects 
for KMM analysis we would have to extend the radial cut to r = 6.'1, thereby including a 
large number of potentially contaminating objects in the sample. We therefore apply the 
empirically-determined cut at B — R — 1.23 and find a blue fraction of 33% (N = 8). 
Increasing the radial cut by 1-2 radial bins (an additional 0^75 — 1'5) does not significantly 
change this result even though the sam ple increases to N = 39 objects. The ACSVCS study 



NGC 4754's g — z color distribution (jPeng et al.1 l2006al ) only detected bimodality at less 



than the 2a significance level (89%) using iV = 83 GC candidates. They estimate a slightly 
larger blue GC fraction of 50%. Although we find no strong statistical evidence for color 
bimodality, we will adopt a blue GC fraction of ~ 40% (average of the empirical split and 
ACSVCS results) in our discussion of the mass-normalized number of blue GCs (T blue ) in 
Section 15.31 Examining the B — R colors for the 90% sample with various radial cuts, the 
lack of sufficient numbers of GC candidates precludes definitive conclusions about a color 
gradient in the NGC 4754 GC system. In addition, we find no evidence for a FF population 
similar to that in NGC 4762. 

NGC 3384. - For NGC 3384, the number of GC candidates in the 90% complete, 
radial cut sample was insufficient to reliably run the KMM tests (N = 13). With no radial 
cut the 90% complete sample only consisted of N = 45 objects, which is still below the 
recommended number of objects (N = 50). We therefore used the empirically-determined 
split at B — R = 1.23 to estimate the blue GC fraction. With no radial cut, the 90% complete 
sample (N = 45) has 30 objects with B — R < 1.23, yielding a blue GC fraction of 67%. 
Relaxing the radial cut imposed at the radial extent of the GC system (r = 5f6; outer edge 
of bin which defines the radial extent) by an additional two radial bins (to r = 7'6) still 
yields only iV = 18 objects, of which 55% (N = 10) are bluer than B — R = 1.23. As our 
final blue GC fraction estimate for NGC 3384, we adopt the average of the estimates from 
the full 90% sample (with no radial cut) and the 90% sample with a radial cut at r = 7'.6, 
finding a value of ~ 60%. Examining the B — R colors for the 90% sample with various 
radial cuts, the lack of sufficient numbers of GC candidates precludes definitive conclusions 
about a color gradient in the NGC 3384 GC system. 
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5.3. Total Numbers of GCs and Global Specific Frequencies 



To estimate the total number of GCs (Nqc) for each galaxy, we integrate the best- 
fit radial surface density profiles as determined from the profile fitting in Section 15.11 For 
each galaxy, we integrate the profile from the inner edge of the innermost bin to the outer 
edge of the bin which defines the radial extent (the bin where the surface density first 
becomes consistent zero within the errors and stays consistent with zero). Because we have 
no information regarding the surface density inwards of the innermost radial bin (due to the 
high surface brightness of the galaxy), we estimate the inner region contribution to iVcc as 
follows. We consider two possible scenarios: that the best-fit profile continues inward to the 
galaxy center or that the profile could remain flat (constant surface density) in the inner 
region. We then take the average of these two estimates as the number of GCs contributed by 
the innermost, unmeasured region. The uncertainty on iVcc includes the following sources 
of error: (1) Poisson errors on the number of GCs and contaminating objects; (2) fractional 
coverage of the GCLF, based on variations due to choices in bin size, dispersions, and bin 
centers; (3) the number of GCs assumed for the central region of the galaxy. For NGC 5866 
and NGC 3384, we include an additional error term due to the larger uncertainty in the 
radial extent of the GC system in these galaxies. 

We calculate the luminosity- and mass-no rmalized specific frequ encies using our esti- 
mates of N GC , mass-to-light ratios (M/L) v from lZepf fc Ashman! (119931 ). and galaxy absolute 
magnitudes from the RC3 (see Tabled]). The luminosity-normalized specific frequency Sn 
i s defined as the total number of GCs normalized by the V band luminosity of the galaxy 
(IHarris & van den Berghlll98ll ). or, 



Sn EE iY GC 10 a4(M 



-15) 



Normalizing the total number of GC s by the stellar mass of the galaxy Mq defines the 
mass- normalized specific frequency T (IZepf fc Ashman! Il993l ): 



T 



GC 



M G /1O 9 M 



(2) 



When calculating values of T, we assume the mass -to- light ratios (M/L) v of 7.6 for lentic- 
ulars and 10.0 for ellipticals (IZepf fc Ashman! 1 19931 ) . In addition to the uncertainty in iVcc 
the errors in Sn and T include an uncertainty in the galaxy magnitude. We adopted a 
galaxy mag nitude uncertainty of 3 time s the standard error in the RC3 magnitude (typically 
~ 0.2 mag; Ide Vaucouleurs et al.lll99ll ) and the uncertainties in all contributing quantities 
were added in quadrature to give the total error in Sn and T. For all galaxies, we find 
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that the uncertainty in the galaxy magnitude is the dominant source of error in the specific 
frequencies. Lastly, we also calculate the mass- normalized number of blue GCs (Tbi ue ), a 
quantity particularly useful in understanding the fi rst generation of GC formation in giant 



galaxies (IRhode et al.ll2005l ; iBrodie &: Straderll2006l ). For each galaxy, we adopt the estimate 
of the fraction of blue GCs determined in Section 15721 for our calculation of Tbiue- 

We discuss the results for the individual galaxies below, comparing our results to pre- 
vious studies where appropriate. Our final values of N GC , Sn, T, and T blue are listed in 
Table [TTJ 

NGC 5866. - For NGC 5866, we integrated the best-fit power law profile from r = 0'.17 
to r = 1V.0 resulting in 268 GCs. Inward of r = Of 17, we estimate an additional 71 GCs 
would contribute to N GC based on the average of a flat profile and the continuation of the 
power law to r = 0. Combining these results yields a final estimate of the total number 
of GCs of iVcc = 340 ± 80. For comparison, CBR07 estimated a total number of GCs of 
~ 400 from their HST study, in good agreement with our results. Using the total absolute 
U-band magnitude and galaxy (stellar) mass from Table [TJ we find specific frequencies of 
Sn = 1-3 ± 0.4 and T = 1.9 ± 0.5. Using the blue GC fraction of 0.74 determined in 
Section [5721 we find a mass-normalized number of blue GCs of Tbi ue = 1.4 ±0.4. Considering 
the various sources of error Nqq, we find that the number of GCs assumed for the inner 
region of the profile is the dominant source of uncertainty in Nqc for NGC 5866. 

NGC 5813. - For NGC 5813, we integrated the best-fit de Vaucouleurs profile from 
r = 0^44 to r = 13M4 resulting in 2743 GCs. Assuming the average of a flat radial profile 
and a continuing de Vaucouleurs profile to r = 0, we estimate an additional 146 GCs would 
contribute to Nqq from the inner, unmeasurable region of the GC system. Combining 
these results gives a total number of GCs Nqq — 2900 ± 400. Because the formal error 
analysis resulted in a fractional uncertainty (t Ngc /Nqq < 10% (due to the relatively low 
uncertainty in the GCLF coverage), we ultimately adopted an error on Nqq commensurate 
with the typical fractional uncertainty from our wide-field imaging survey (~ 10% — 30%). 
Using the total absolute l/-band magnitude and galaxy (stellar) mass from Table [fl we find 
specific frequencies of Sn = 3.6 ± 0.8 and T = 4.2 ± 0.9. Using the blue GC fraction of 0.68 
de termined in Section 15.21 we find a mass-normalized number of blue GCs of Tbi ue = 2.8 ±0.6 



Harris fc van den Berghl (1198ll ) found Nqq = 2400 ± 700, and although th eir result is 



consis tent with ours, their uncertainties are nearly two times as large as ours. iHopp et al. 



(119951 ) estimated Nqq = 1650 ±400, approximately 3a smaller than our estimate. They find 
that the radial surface density of their GC candidates extends to only a projected radius of 
3'.2 (compared to our estimate of r ~ 13'), likely contributing to their lower estimate of Nqq. 



NGC 4^62. - For NGC 4762, we integrated the best-fit de Vaucouleurs profile from 
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r = O'.l to r = 5'.5 resulting in 261 GCs. Inwards of r = Oil, we estimate that an additional 
9 GCs would contribute to the total based on the average of a flat inner profile and the 
continuation of the de Vaucouleurs profile inward to r = 0. We therefore find a total number 
of GCs of iVcc = 270 ± 30, where the dominant source of error is the uncertainty in the 
GCLF correction. We find specific frequencies of Sn = 0.9 ± 0.2 and T = 1.4 ± 0.3 using 
our estimate of Ago the galaxy absolute V^-band magnitude, and the galaxy (stellar) mass 
from Table [TJ Using the blue GC fraction of 0.40 determined in Section 15.21 we find a 



mass- normalized number of blue GCs of Tbi uc = 0.6 ± 0.1. For comparison, iPeng et al. 



( 120081 ) estimate a total nu mber of GC s in NG C 4762 of 211 ± 34 from their analysis of the 



ACSVCS data. Although IPeng et al.l (120081 ) do include an extrapolation of their best-fit 



Sersic radial surface density profile to account for the limited ACS/WFC field of view, it 
is unclear to what radius they integrated the profile to derive Ago Examining our radial 
profile integration results, we find that the surface density profile outside a projected radius 
of ~ 2f7 (the limit of the ACS/WFC) contributes ~ 50 GCs to A GC - nearly 20% of the 
total. Not surprisingly then, differences and/or uncertainties in the assumed radial extent 
and surface density profile shapes can contribute significantly to systematic uncertainties in 
the total number of GCs. 

NGC 4^54- - For NGC 4754, we integrated the best-fit de Vaucouleurs profile from 
r = 0'.09 to r = 3'5 resulting in 113 GCs. Inwards of r = 0'.09, we estimate that an 
additional 2 GCs would contribute to the total based on the average of flat inner profile 
and the continuation of the de Vaucouleurs profile inward to r = 0. We therefore find a 
total number of GCs of Nqq = 115 ± 15, where the dominant source of uncertainty is the 
Poisson errors on the number of GCs and contaminating objects. We find specific frequencies 
of Sn = 0.6 ±0.1 and T = 0.9 ± 0.2 using our estimate of Aqc and the galaxy absolute 
magnitude and (stellar) mass from Table [TJ Using the blue GC fraction of 0.40 determined 
in Section 15.21 we find a mass-norm alized blue GC fra ction of Tbi ue = 0.4 ± 0.1. Based on 



their ACSVCS study of NGC 4754, IPeng et al.l (120081 ) estimate a total number of GCs in 



NGC 4754 of 103 ±17 (based on an extrapolated Sersic fit), in good agreement with our 
results. 

NGC 3384. - For NGC 3384, we integrated the best-fit power law profile from r = 0'.59 
to r = 5f6 resulting in 106 GCs. Inwards of r = 0'.59, we estimate than an additional 16 
GCs would contribute to the total based on the average of the flat inner profile and the 
continuation of the best-fit power law profile inward to r = 0. We thus find a total number 
of GCs Nqq = 120 ± 30. We found that all terms in the error analysis contributed nearly 
equally; no source of uncertainty dominates the error in N GC . We find specific frequencies 
of Sn = 0.8 ± 0.2 and T — 1.2 ± 0.4 using our estimate of Aqc and the galaxy absolute 
magnitude and (stellar) mass from Table [TJ Using the blue GC fraction of 0.60 determined in 
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Section 15.21 we find a mass-nor malized blue GC fraction of Tbi uo = 0.7±0.2. For comparison, 
Harris & van den Berghl (Il98ll ) estimated Nqq = 110 ±60 for NGC 3384 from their analysis 



of the N3379 field. Although their photographic study is in good agreement with our result, 
our three-filter CCD observations and analysis methods have reduced the uncertainty in Nqq 
by a factor of two. 



5.3.1. Comparison of Results to Other Giant Galaxies 

A primary goal of our extragalactic GC system studies is to better quantify how the 
global properties of GC systems vary with host galaxy properties such as luminosity, envi- 
ronment, and morphology. Our analysis of the five early- type galaxies in this study adds one 
elliptical and four lenticular galaxies to our larger wide-field imaging survey. This doubles 
the number of lenticular galaxies in our survey sample and results in a comparable num- 
bers of ellipiticals and SOs. In combination with the spiral galaxies, our current sample of 
giant galaxies with well-measured total numbers of GCs and specific frequencies consists of 
26 galaxies: 11 spirals, 8 lenticulars, and 7 ellipticals. Twenty-one of these are from our 
wide-field imaging survey and five galaxies are from the literature, including the Milky Way 
and M31. The results for 22 of these galaxies are summarized in Y12 and this work adds an 
additional four galaxies. (The results for NGC 3384 were included in the analysis by Y12.) 
Literature references and values for N GC , Sn, T, and T blue for the 22 galaxy sample can be 
found in Table 10 of Y12. Although future study will perform a more detailed analysis of 
the full survey sample (Hargis & Rhode 2013, in preparation), here we compare values of 
Sjv and T by host galaxy morphological type. 

To investigate the variation in specific frequencies with galaxy morphology, we calculated 
the weighted mean values of Sn and T for the sample of 26 elliptical, spiral, and lenticular 
galaxies. The sample of spiral galaxies is identical to that presented in Y12; for 11 galaxies 
(including the Milky Way and M31) they find weighted mean specific frequencies of Sn = 
0.6 ± 0.1 and T = 1.2 ± 0.1. For the sample of seven elliptical galaxies we find weighted 
means of Sn = 1.8 ± 0.1 and T = 2.1 ± 0.2, consistent with the combined analysis of 11 
E/S0 galaxies by Y12. For the sample of eight lenticular galaxies we find weighted means 
of Sn = 0.9 ± 0.1 and T = 1.5 ± 0.1. Consistent with previous results on giant galaxy 
GC systems, we find that ellipticals have larger luminosity- and mass-normalized specific 
frequencies than spirals. Comparing the lenticulars to spirals and ellipticals, we find that 
on average, SOs have GC specific frequencies more comparable to spirals than to ellipticals. 
The difference in mean values of Sn and T between spirals and lenticulars (~ 0.3) is only 
marginally larger than the formal uncertainties in the weighted means; the mean values are 
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consistent within 3a. For the galaxies in this work, we find that the four lenticulars occupy 
a range of specific frequencies consistent with both the weighted means of spirals and SOs. 
Quantifying the degree to which the GC systems of SOs and spirals are statistically different 
(or similar) will require a larger sample of galaxies with well-measured global properties. 

6. Summary and Conclusions 

In this study we investigated the global properties of the GC systems of five moderate- 
luminosity {My ~ —21 to —22), early- type giant galaxies from our wide- field imaging survey 
program. We used KPNO 4-m MOSAIC BVR imaging to trace the GC systems to large 
projected radii. The final results of our study are given in Table ITT1 Here we summarize our 
main results: 

1. We constructed radial surface density profiles for the GC systems for the target galaxies 
using the MOSAIC imaging data and, where available, complementary archival and 
published HST data. While the large spatial coverage of the MOSAIC imager allows us 
to trace the GC population to galactocentric radii of ~ 150 — 200 kpc (for the distances 
of our target galaxies), the higher spatial resolution of HST instruments allows us to 
probe the inner several hundred parsecs. 

2. We integrate the best-fit radial surface density profiles to derive estimates of the total 
number of GCs {Nqq). We find total numbers of GCs ranging from Agc ~ 120 for 
the two lowest luminosity galaxies to Nqq ~ 2900 for the most luminous galaxy in our 
study. We find C-band normalized specific frequencies Sn ranging from 0.6 to 3.6 and 
stellar mass-normalized specific frequencies of T ranging from 0.9 to 4.2. Analyzing a 
larger sample of 26 galaxies with well-measured total numbers and specific frequencies, 
we find that spiral and lenticular galaxies have statistically consistent (weighted) mean 
values of Sn and T. The four SO galaxies in this study have values of Sjy and T 
comparable to the mean values of both spiral and lenticular galaxies. 

3. We studied the color distributions and color gradients of the galaxies' GC systems 
using the MOSAIC imaging data. Only NGC 5813 shows strong evidence (99.99% 
confidence level) of color bimodality. NGC 5866 shows some evidence for color bi- 
modality, although at a lower level of statistical significance (97% confidence level). 
NGC 5813 has sufficient numbers of GC candidates to investigate the spatial distribu- 
tions of the red (metal-rich) and blue (metal-poor) GC subpopulations. The red GC 
subpopulation is more centrally concentrated about the galaxy center than the blue 
GC subpopulation and a statistically significant B — R color gradient is found in the 
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inner 6' (~ 55 kpc) of the GC system. No statistically significant color gradients are 
found in the GC systems of the other target galaxies. 

The primary goal of this work has been to measure the global properties of the GC 
systems for a number of intermediate-luminosity, early-type giant galaxies. While a fair 



number of t 
techniques ( 


re most massive early-type galaxies have been studied using wide-field imaging 


Dirsch et a 


2003: 


Harris et al. 


2004; 


Rhode & Zepi 


2001. 


2004; 


Tamura et al. 


2006; 


Forbes et al. 


2011) 


, studies with large spatial coverage are rare for moderate-luminosity 



galaxies. With the addition of these five galaxies to our larger survey sample, a future study 
(Hargis & Rhode 2013, in preparation) will present the results of a multi-variate analysis of 
the current survey sample of ~ 30 galaxies. 
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Fig. 1. — Digitized Sky Survey image of NGC 5866 showing the area covered by the Mosaic 
imaging analyzed in this study (large box). The HST ACS pointing analyzed in this study 
is also shown (smaller boxes). 




Fig. 2. — Digitized Sky Survey image of NGC 5813 showing the area covered by the Mosaic 
imaging analyzed in this study (large box). 




Fig. 3. — Digitized Sky Survey image of the NGC 4754/4762 field showing the area covered 
by the Mosaic imaging analyzed in this study (large box). The two HST ACS pointings 
analyzed in this study are also shown (smaller boxes). 
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Fig. 4. — Digitized Sky Survey image of the NGC 3384/3379 field showing the area covered 
by the Mosaic imaging analyzed in this study (large box). 
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Fig. 5. — Example point source selection for the detected objects in the Mosaic images. 
Plotted are the full-width at half-maximum (FWHM) of the Gaussian profiles versus instru- 
mental BVR magnitudes of detected sources which were selected as point sources (filled 
circles) and rejected as extended objects (open squares). Extended source cuts similar to 
this example were done for all target galaxies. 
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Fig. 6. — Selection of the Mosaic GC candidates in NGC 5866 based on V magnitudes and 
B — V, V — R colors. The point sources detected in all three filters are shown as open 
squares and the final sample of 290 GC candidates are shown as filled circles. The red box 
denotes the boundary of the color selection (see Section I3T21 . The imposed V magnitude cut 
rejects some point sources in the selection box. Objects outside the selection box included as 
GC candidates results from the consideration of photometric errors in the selection process. 
The lower right corner shows the median B — V and V — R photometric errors for the GC 
candidates. The curves denote the location of galaxies of various morphological types from 
redshift z = to 0.7. The upper left corner shows the reddening vector corresponding to 

Ay = 1. 
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Fig. 7. — Selection of the Mosaic GC candidates in NGC 5813 based on V magnitudes and 
B — V, V — R colors. The point sources detected in all three filters are shown as open 
squares and the final sample of 1300 GC candidates are shown as filled circles. See caption 
for Figure [6] for additional plot details. 
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Fig. 8. — Selection of the Mosaic GC candidates in NGC 4762 based on V magnitudes and 
B — V, V — R colors. The point sources detected in all three filters are shown as open 
squares and the final sample of 481 objects (GC candidates) with V magnitude and BVR 
colors that meet our selection criteria are shown as filled circles. The selection has been done 
for the entire Mosaic image (~ 36' x 36') with the region around NGC 4754 masked out. See 
caption for Figure [6] for additional plot details. 
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Fig. 9. — Selection of the Mosaic GC candidates in NGC 4754 based on V magnitudes and 
B — V, V — R colors. The point sources detected in all three filters are shown as open 
squares and the final sample of 451 objects (GC candidates) with V magnitude and BVR 
colors that meet our selection criteria are shown as filled circles. The selection has been done 
for the entire Mosaic image (~ 36' x 36') with the region around NGC 4762 masked out. See 
caption for Figure [6] for additional plot details. 
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Fig. 10. — Selection of the Mosaic GC candidates in NGC 4754 based on V magnitudes 
and B — V, V — R colors. The point sources detected in all three filters are shown as open 
squares and the final sample of 181 GC candidates are shown as filled circles. See caption 
for Figure [6] for additional plot details. 
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Fig. 11. — Color-color diagram for a Besangon Galactic star count model for the NGC 5866 
Mosaic field. Estimates of the foreground Galactic contamination were made using the 
Besangon stellar population synthesis models (see Section 14.2. 2p . The V magnitude range of 
the model stars was chosen to match the GC candidate selection range for NGC 5866 (see 
Section 13^21) . Open squares show the B — V and V — R colors of all model stars. The scatter 
in the colors is intrinsic to the simulations. Solid circles show the model stars that would 
pass the color selection criteria (defined by the box) for our GC candidates and would be 
considered contaminating objects. The background galaxy curves and reddening vector are 
identical to Figures I&HT01 and are shown for comparison to these diagrams. 
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Fig. 12. — GCLF fits for the target galaxies. The observed GCLF is shown as the shaded 
histogram and the completeness-corrected GCLF (used in the fits) is shown as the solid line 
histogram. The best-fit theoretical Gaussian functions are shown for three values of the 
dispersion: a = 1.2 (dotted line), 1.3 (short dashed line), and 1.4 (long dash). The arrow 
denotes the magnitude where the convolved BVR completeness is 50%. 
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Fig. 13. — Radial profile of the NGC 5866 GC system. The profile has been corrected 
for missing area, contamination, and magnitude incompleteness (see Section IBTTj) . The top 
panel shows the corrected surface density as a function of projected radial distance from the 
galaxy center (r = 0). The bottom panel shows the corrected surface density as a function 
of projected radius on a logarithmic scale. The filled circles are the Mosaic data points 
from this study and the fi lled triangles are our analysis of the HST ACS GC candidates 
from ICantiello et al.l (120071 ) . Both panels show the best fit power-law and r 1 / 4 profile fits 
(described in Section 15. ip as the solid line and dashed line, respectively. 
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Fig. 14. — Radial profile of the NGC 5813 GC system plotted in the same manner 
Figure [13j See caption for Figure [13] for additional plot details. 
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Fig. 15. — Radial profile of the NGC 4762 GC system plotted in the same manner as 
Figure [13j See caption for Figure [13] for additional plot details. 
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Fig. 16. — Radial profile of the NGC 4754 GC system plotted in the same manner as 
Figure [13j See caption for Figure [13] for additional plot details. 
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Fig. 17. — Radial profile of the NGC 3384 GC system plotted in the same manner 
Figure [13j See caption for Figure [13] for additional plot details. 
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Fig. 18. — Color Distribution and KMM mixture modeling results for NGC 5866. The shaded 
histogram shows the B — R color distribution for the radial cut (r < 11'), 90% magnitude 
complete sub-sample of 51 Mosaic GC candidates for NGC 5866 (see Section 15. 2j) . The 
KMM mixture modeling results are shown as the solid and dashed curves. The solid curve 
denotes the blue (metal-poor) population and the dashed curve denotes the red (metal-poor) 
population. 
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Fig. 19. — B — R color versus projected radial distance of the 90% sample of GC candidates 
for NGC 5866, NGC 5813, and NGC 4762. GC candidates in the 90% sample are shown 
as open circles. For NGC 5866 and NGC 4762, we performed linear fits to the data using 
only the GC candidates within the radial extent of the GC system and find no evidence of 
statistically signficant color gradients. For NGC 5813, we performed linear fits over the full 
radial extent (r < 13f4) and in the inner regions of the GC system. We find evidence of 
a statistically significant color gradient in the inner 6' (~ 55 kpc) of the GC system (slope 
= A(B — R)/Ar = —0.018 ± 0.006). For each galaxy we show the linear fit as a red solid 
line and the best-fit slope and la uncertainty. NGC 4754 and NGC 3384 had insufficient 
numbers of GC candidates in the 90% sample (within the radial extent of the GC system) 
so no color gradient analysis was performed. 
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Fig. 20. — Color Distribution and KMM mixture modeling results for NGC 5813. The shaded 
histograms (both panels) show the B — R color distribution for the radial cut (r < 13!44), 
90% magnitude complete sub-sample of 528 Mosaic GC candidates (see Section l5\2l . The 
KMM mixture modeling results for the homoscedastic (same dispersions) case in panel (a) 
and the results for the heteroscedastic case (differing dispersions) case in panel (b). In both 
panels the solid curves denote the blue (metal-poor) population and the dashed curves denote 
the red (metal-poor) population. 
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Fig. 21. — Spatial distribution of the red and blue GC subpopulations in NGC 5813. The 
upper panels shows the spatial positions of the red (metal-rich) and blue (metal-poor) GC 
candidates on the Mosaic frame as blue filled circles (upper left panel) and red filled circles 
(upper right panel). The center of NGC 5813 is shown as the black cross. The radial extent 
of the full GC system (red+blue subpopulations; see Figure [H} is shown as the circle of 
radius r = 12'.9 (120 kpc). The lower panel shows the radial surface density profiles of the 
two subpopulations. We show only the inner 15' of the profile; the complete radial profile 
for the full (blue+red) GC population is shown in Figure HU The red GC subpopulation is 
shown as the red open circles and the blue subpopulation is shown as the blue filled circles. 
The spatial distributions of the subpopulations are discussed in Section 15.21 



Table 1. Basic Properties of Target Early- Type Galaxies 



Name (VCC ID) 


Type 


Ml 


m — M 


Distance (Mpc) 


vs 


Ay 


"holio (km s 1 ) 


log M/Mq 


Environment 


NGC 5866 


SAO 


-21.1 


30.93 ±0.12 


15.3 


9.85 1 


0.044 


672 


11.2 


NGC 5866 Group 


NGC 5813 


El 


-22.3 


32.54 ±0.18 


32.2 


10.27 1 


0.189 


1972 


11.8 


NGC 5846 Group 


NGC 4762 (VCC 2095) 


SBO 


-21.2 


31.35 ±0.15 2 


18.6 


10.16 


0.086 


984 


11.3 


Virgo Cluster 


NGC 4754 (VCC 2092) 


SBO 


-20.7 


31.13 ±0.14 


16.8 


10.43 


0.086 


1347 


11.1 


Virgo Cluster 


NGC 3384 


SBO 


-20.5 


30.32 ±0.14 


11.6 


9.84 


0.081 


704 


11.0 


Leo I Group 



Note. — Galaxy morpho logical types are from Ide Vaucouleurs e^_al ] lll99ll) . MZ are computed by combining m — M with V^?. V£ are from 
deVaucouleursetalJ ill 9911) except where noted, m — M are froni iTMrrveTaLl l l200ll ). except where noted. Galactic extinction values Ay are from 
SchIegel*et*^lTTl998l ) . Heliocentric radial velocities ^ helio are from the NASA /IPAC Extragalactic Database (NED). Galaxy stellar masses log M/Mq 
are computed from My and assumed (M/L)y from lZepf fc Ashman! il993T l of 7.6 for lenticulars and 10.0 for ellipticals. 

1 Computed from V~t from lde"V aucoulcur s et al. (1991) and the adopted Ay. 

2 Computed by combining the 3K CMB recession velocity of 1302 ± 25 km s" 1 (NED) with H = 70 ± 5 km s^ 1 Mpc" 1 , 
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Table 2. Mosaic Observations of Target Galaxies. 



Name 


Date 


Exposure Times (s) 










B 


V 




R 


NGC 5866 


May 2010 


4 x 800; 1 x 600 


5 x 700 


3 


x 600 


NGC 5813 


May 2010 


4 x 1260 


5 x 1100 


4 


x 1080 


NGC 4762/4754 


May 2010 


5 x 1260 


4 x 1200 


3 


x 1080 


NGC 4216 


May 2010 


5 x 1260 


5 x 1200 


5 


x 1080 


NGC 3384 


March 1999 


5 x 660 


5 x 540 


5 


x 480 



Table 3. 50% and 90% Completeness Limits for Mosaic Data 



Image NGC 5866 NGC 5813 NGC 4762/4754 NGC 3384 



50% 90% 50% 90% 50% 90% 50% 90% 



B 23.80 23.46 25.13 24.80 24.63 24.13 23.56 22.60 

V 24.03 23.63 24.23 23.85 23.97 23.62 24.09 23.19 

R 23.46 23.03 24.07 23.48 23.17 22.82 23.40 22.84 
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Table 4. HST Observations Analyzed for this Study 



Proposal ID 


PI 


Instrument Exposure Time (s) 


Filter 


NGC 5866 


10705 


Noll 


ACS/WFC 


3900 


F435W 


10705 


Noll 


ACS/WFC 


2800 


F555W 


10705 


Noll 


ACS/WFC 


2200 


F625W 


978 


Ho 


ACS/WFC 


120 


F814W 


NGC 5813 


5454 


Franx 


WFPC2 


1000 


F555W 


5454 


Franx 


WFPC2 


460 


F814W 


6357 


Jaffe 


WFPC2 


1000 


F702W 


NGC 4762 


9401 


Cote 


ACS/WFC 


750 


F475W 


9401 


Cote 


ACS/WFC 


1210 


F850LP 


NGC 4754 


9401 


Cote 


ACS/WFC 


750 


F475W 


9401 


Cote 


ACS/WFC 


1210 


F850LP 



Note. — HST observations of NGC 3384 were analyzed in RZ04. 
See their Table 5 for details. 
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Table 5. Corrected Radial Surface Density Profile for GC Candidates in NGC 5866 



Projected Radius (arcmin) 


Surface Density (arcmin 2 ] 


Fractional Coverage 


Source 


0.4 


32.46 ± 7.26 


0.89 


HST 


0.8 


16.55 ±3.38 


1.00 


HST 


1.2 


6.81 ± 1.76 


0.99 


HST 


1.6 


3.90 ±1.30 


0.77 


HST 


2.0 


1.88 ±1.33 


0.28 


HST 


2.6 


1.42 ±0.41 


0.82 


Mosaic 


4.9 


0.66 ±0.20 


1.00 


Mosaic 


7.3 


0.23 ±0.14 


0.94 


Mosaic 


9.8 


0.13 ±0.12 


0.78 


Mosaic 


12.3 


0.08 ±0.11 


0.72 


Mosaic 


14.8 


0.06 ±0.11 


0.65 


Mosaic 


17.2 


0.09 ±0.11 


0.50 


Mosaic 


19.7 


-0.10 ±0.11 


0.33 


Mosaic 



Note. — Negative surface densities can occur due to the application of the contamination 
correction. 
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Table 6. Corrected Radial Surface Density Profile for GC Candidates in NGC 5813 



Projected Radius (arcmin) 


Surface Density (arcmin 2 " 


Fractional Coverage 


Source 


1.1 


77.27 ±8.82 


0.90 


Mosaic 


2.0 


35.88 ±4.18 


0.94 


Mosaic 


3.0 


18.51 ±2.66 


0.86 


Mosaic 


4.0 


10.24 ± 1.73 


0.96 


Mosaic 


5.0 


6.93 ±1.35 


0.96 


Mosaic 


6.0 


5.38 ±1.14 


0.95 


Mosaic 


7.0 


3.78 ±0.96 


0.96 


Mosaic 


8.0 


2.44 ±0.82 


0.96 


Mosaic 


8.9 


1.72 ±0.74 


0.92 


Mosaic 


10.0 


0.69 ±0.63 


0.95 


Mosaic 


10.9 


1.36 ±0.64 


0.97 


Mosaic 


11.9 


0.81 ±0.59 


0.94 


Mosaic 


12.9 


0.27 ±0.53 


0.96 


Mosaic 


13.9 


1.02 ±0.58 


0.87 


Mosaic 


14.9 


-0.09 ±0.52 


0.79 


Mosaic 


15.9 


0.47 ±0.54 


0.78 


Mosaic 


16.9 


0.22 ± 0.52 


0.76 


Mosaic 


17.9 


0.51 ±0.56 


0.64 


Mosaic 


18.9 


-0.86 ±0.51 


0.50 


Mosaic 


19.9 


0.27 ±0.64 


0.41 


Mosaic 


20.9 


-0.18 ±0.65 


0.34 


Mosaic 


21.9 


-0.44 ±0.65 


0.29 


Mosaic 


22.9 


0.08 ±0.82 


0.20 


Mosaic 


23.9 


0.35 ±0.98 


0.14 


Mosaic 


24.9 


0.32 ±1.12 


0.10 


Mosaic 


25.9 


-0.69 ±1.17 


0.07 


Mosaic 



Note. — Negative surface densities can occur due to the application of the contamination 
correction. 
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Table 7. Corrected Radial Surface Density Profile for GC Candidates in NGC 4762 



Projected Radius (arcmin) 


Surface Density (arcmin 2 ) 


Fractional Coverage 


Source 


0.5 


42.62 ± 7.03 


0.97 


HST 


1.0 


15.44 ±2.69 


0.99 


HST 


1.3 


9.08 ±2.21 


0.77 


Mosaic 


1.6 


5.68 ±1.55 


0.74 


HST 


2.1 


3.05 ±1.88 


0.21 


HST 


2.4 


1.55 ±0.79 


0.79 


Mosaic 


3.7 


0.78 ±0.52 


0.88 


Mosaic 


4.9 


0.15 ±0.37 


0.94 


Mosaic 


6.2 


0.17 ±0.33 


0.97 


Mosaic 


7.4 


-0.12 ±0.28 


0.92 


Mosaic 


8.6 


0.62 ±0.33 


0.88 


Mosaic 


9.9 


0.15 ±0.28 


0.87 


Mosaic 



Note. — Negative surface densities can occur due to the application of the contamination 
correction. 
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Table 8. Corrected Radial Surface Density Profile for GC Candidates in NGC 4754 



Projected Radius (arcmin) 


Surface Density (arcmin 2 ) 


Fractional Coverage 


Source 


0.2 


53.00 ± 18.93 


1.00 


HST 


0.5 


16.22 ±3.85 


0.97 


HST 


0.9 


5.20 ±2.39 


0.81 


Mosaic 


1.4 


5.05 ±0.86 


0.82 


HST 


1.6 


1.70 ±1.07 


0.95 


Mosaic 


2.3 


0.79 ±0.73 


0.94 


Mosaic 


3.1 


0.33 ±0.56 


0.94 


Mosaic 


3.8 


0.02 ±0.45 


0.96 


Mosaic 


4.6 


0.49 ±0.48 


0.97 


Mosaic 


5.3 


-0.01 ±0.39 


0.91 


Mosaic 


6.1 


0.28 ±0.42 


0.85 


Mosaic 


6.8 


0.36 ±0.41 


0.83 


Mosaic 


7.6 


0.85 ±0.45 


0.81 


Mosaic 


8.3 


0.29 ±0.38 


0.76 


Mosaic 



Note. — Negative surface densities can occur due to the application of the contamination 
correction. 
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Table 9. Corrected Radial Surface Density Profile for GC Candidates in NGC 3384 



Projected Radius (arcmin) 


Surface Density (arcmin 2 ) 


Fractional Coverage 


Source 


1.2 


3.48 ± 1.24 


0.85 


Mosaic 


2.1 


1.70 ±0.61 


0.98 


Mosaic 


3.1 


1.36 ±0.52 


0.79 


Mosaic 


4.1 


0.34 ±0.31 


0.70 


Mosaic 


5.1 


0.00 ±0.22 


0.68 


Mosaic 


6.1 


0.03 ±0.21 


0.66 


Mosaic 


7.1 


0.11 ±0.20 


0.68 


Mosaic 


8.1 


0.23 ±0.21 


0.69 


Mosaic 


9.1 


0.03 ±0.17 


0.71 


Mosaic 


10.1 


0.17 ±0.18 


0.72 


Mosaic 


11.1 


-0.02 ±0.14 


0.73 


Mosaic 


12.1 


0.17 ±0.16 


0.77 


Mosaic 


13.1 


-0.02 ±0.13 


0.76 


Mosaic 


14.1 


-0.08 ±0.12 


0.75 


Mosaic 


15.1 


0.25 ±0.16 


0.72 


Mosaic 


16.1 


-0.13 ±0.11 


0.61 


Mosaic 


17.1 


0.11 ±0.15 


0.56 


Mosaic 


18.1 


-0.11 ±0.12 


0.51 


Mosaic 



Note. — Negative surface densities can occur due to the application of the contamination 
correction. 



Table 10. GC Radial Surface Density Profile Fit Parameters 



Power Law de Vaucouleurs Law 



Galaxy 


a 


ai 


X 2 /" 


a 


ai 


X 2 /" 


NGC 5866 


0.96 ±0.04 


-1.76 ±0.09 


0.47 


3.52 ±0.16 


-2.52 ±0.14 


0.70 


NGC 5813 


2.02 ±0.04 


-1.74 ±0.07 


0.86 


4.18 ±0.12 


-2.23 ±0.09 


0.33 


NGC 4762 


1.10 ±0.04 


-1.86 ±0.14 


1.02 


4.23 ±0.23 


-3.09 ±0.23 


0.52 


NGC 4754 


0.78 ±0.05 


-1.41 ±0.13 


1.36 


3.44 ±0.24 


-2.61 ±0.24 


0.97 


NGC 3384 


0.65 ±0.11 


-1.29 ±0.17 


0.88 


2.06 ±0.28 


-1.51 ±0.20 


1.02 



Table 11. Total Numbers of GCs and Specific Frequencies of GC Systems for Target Galaxies 



Galaxy 


N GC 




T 


Thine 


/blue 


Color Bimodality 


Color Gradient 


Radial Extent 














(Yes/No) 


(Yes/No) 


(kpc) 


NGC 5866 


340 ± 80 


1.3 ±0.4 


1.9 ±0.5 


1.4 ±0.4 


0.74 


Y 


N 


44 ± 11 


NGC 5813 


2900 ± 400 


3.6 ±0.8 


4.2 ±0.9 


2.8 ±0.6 


0.68 


Y 


Y 


120 ± 14 


NGC 4762 


270 ± 30 


0.9 ±0.2 


1.4 ±0.3 


0.6 ±0.1 


0.40 


N 


N 


27 ±7 


NGC 4754 


115 ± 15 


0.6 ±0.1 


0.9 ±0.2 


0.4 ±0.1 


0.40 


N 


N 


15 ±4 


NGC 3384 


120 ± 30 


0.8 ±0.2 


1.2 ±0.4 


0.7 ±0.2 


0.60 


N 


N 


17 ±4 



